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The growth of a rough interface through a random media is 
modelled by a continuous stochastic equation with a quenched 
noise. By use of the Novikov theorem we can transform the 
dependence of the noise on the interface height into an effec- 
tive temporal correlation for different regimes of the evolution 
of the interface. The exponents characterizing the roughness 
of the interface can thus be computed by simple scaling ar- 
guments showing a good agreement with recent experiments 
and numerical simulations. 
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The growth of rough interfaces through a random 
medium is an interesting pattern-formation phenomena 
that has received a lot of attention in recent years [p]-p2[ . 
In particular, experimental studies of the displacement of 
a nonwetting fluid by a more viscous wetting one through 
a porous media have shown that the interface separating 
both fluids is rough and the random pinning of the dis- 
placing fluid on the wider porous of the medium strongly 
influences the growth |^,^. The random media is frozen 
in time and a natural way to model its effects is by in- 
troducing a quenched disorder in a continuous equation, 
as has been done by Kessler, Levine and Tu (KLT) 
Parisi [Q , and others • This model is different from 
the one that describes the growth of solids from a va- 
por or from a stochastic ballistic deposition, in which the 
stochastic field changes both in space and time. These 
last phenomena are usually described by equations of the 
type introduced by Edwards and Wilkinson (EW) fl^ j 
and by Kardar, Parisi and Zhang (KPZ) Q. The ef- 
fects of a quenched noise in a stochastic model are much 
less known and are the subject of this Letter. 

We study the model introduced by KLT |^ and Parisi 
0. Starting from an initially flat interface above a d- 
dimensional substrate at time i = 0, a rough interface is 
described by its height h{x, t) at position x and time t, 
which obeys the following stochastic equation: 
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h{x, t) = V^h{x, t) + F + r]{x, h), 



(1) 



where the diffusive term models the surface tension ef- 
fects, F is the pushing force and the noise 77 is Gaussian, 



< V{x,y)v{x',y') >= -S'^ix-x'jSiy-y'), 

with an intensity 9. This equation also models the growth 
of domain walls in a random-field Ising ferromagnet [p^ . 

Analytical results for the exponents associated with 
the model described by Eq. (|^) arc only known in the 
limit of strong pushing, where the interface moves very 
fast (for d — 1 the roughening exponent a = 1/2 |^ and 
the time exponent /3=l/4 0). In this limit, the model 
reduces to the case of a noise with a delta correlation in 
time and the exponents correspond to the EW model, for 
which the scaling relation between exponents is given by 
a/f3 = 2. On the other hand for small pushing {F < Fc) 
the interface becomes pinned. 

KLT and Parisi have performed numerical sim- 
ulations of the model given by Eq. (|l|) for different values 
of the pushing force. Far from Fc, Parisi Q has ob- 
tained that, after a transient, there is a temporal regime 
with anomalous exponents followed by a crossover to an 
asymptotic regime that corresponds to the EW model. 
As one approaches the pinning transition, the crossover 
appears at later times, but the long-time behavior is de- 
scribed in any case by the EW model. Only at the pin- 
ning transition one gets a different exponent, /3 — 3/A 
0. In the intermediate temporal regime, the numerical 
results of KLT |] give a = 0.73. These results should be 
compared with the experimental results of Rubio et al. 
§ (a = 0.73 ±0.03) and Horvath et al. § (/? = 0.65 and 
a ~ 0.81). 

The purpose of this Letter is to obtain analytical ex- 
pressions for the exponents characterizing the interface 
growth and the associated scaling relations in regimes 
where scaling properties are obeyed. To do so, we start 
with Eq. (P and we derive an effective model valid for 
small intensities of the noise, 9. As we will see below, 
the resulting effective model is given by the following 
stochastic linear differential equation 



dy 
dt 



(2) 



where y{x,t) = h{x,t) — h{t) measures the fluctuations 
of the height interface h{x, t) around its averaged value 
h(t). i,{x,t) is the effective noise with zero mean value 
and correlation given by: 
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< a^,mx',t') -5''{s-x')w{t - 1'). (3) 

The noise is delta-correlated in the substrate position 
and the dependence on h is substituted by a temporal 
correlation W{t — t'). The explicit form of W{t — t') will 
be discussed in the following paragraph. 

From Eqs. it is possible to derive expressions for 

the width of the interface, a —< and the correla- 

tion lenght Ic in terms of W{t — t'). Therefore, the scaling 
properties of W{t — t') will determine the roughening ex- 
ponents and the associated scaling relations. 

Derivation of the effective model.- In order to get an 
effective model, we use functional methods (a generaliza- 
tion [|l6| of the Novikov theorem ||l^) to determine the 
mean value of rj: 



< T]{x, h{x, t)) >= X / G(0, t - t')dt' X 
2 Jo 

X JdyJ dy'S'iy - y')Ps,t;S,t' {y, y') + 0{e^) (4) 
and also its correlation 

Q 

< T]{x,h{x,t))r]{x',h{x',t')) >= -S'^{x-x')x 



X / dyPs^f^s,t'{y,y) + 0{9^) 



(5) 



where G{x,t) ^ t '^^^exp{^-) is the diffusive propaga- 
tor and 

Ps.t;S',t'iy,y') =< 5{y-h{x,t))5{y' -h[x',t')) > (6) 

is the joint probability for the interface to have heights 
y and y' at x,t and x',t' respectively. We note that 
Eqs.(^ and ^ are valid for small values of the intensity 
of the noise because we have neglected terms of order 
9"^ involving the joint probability of four points. By as- 
suming homogeneity both in space and time and defining 
W{u, t — t') as the probability of having an increment u 
in a time t - t', W{u, t-t') ^ Px.us^t' (u + h, hVPs^t' (h), 
we obtain the effective model given by Eqs.(g,^ where 
^{x,t) = r]{x,h{x,t))~ < r]{x,h{x,t)) > is the effec- 
tive noise with zero mean value and correlation given 
by Eq.(^) up to the lowest order in 9, and W{t — t') = 
W{0, t — t') is the probability for the interface to remain 
at the same point after a time interval {t — t'). Now, 
since the model is linear, the kinetics of the interface is 
completely characterized by this probability. 

Exponents and scaling relations.- First, we will dis- 
cuss the different relevant length scales associated with 
the interfacial growth. The interplay between these scales 
will determine the different scaling regimes. One scale is 
given by the horizontal correlation length, Zc, which for 
a very general form of W{t — t'), has diffusive behavior 



with time, Ic t^^"^. This result is remarkable because it 
indicates a very universal behavior for Ic in the complete 
temporal regime independent from the properties of the 
disorder. It only depends on the diffusive character of 
Eq.(||) and on the delta spatial correlation of the noise. 
Another relevant scale is associated to the characteristic 
length, Lc, of the clusters of points of the interface that 
become pinned by the effect of the quenched disorder. 
This cluster length depends on the pushing force. At the 
pinning transition, with a pushing force Fc, the entire in- 
terface becomes pinned and the cluster length is Lc L. 
By increasing F, Lc is reduced and it becomes zero in 
the limit of strong pushing, where the interface moves 
very fast and the model reduces to EW model. In this 
limit, W{t — t') becomes a delta function. In accordance 
with these results, we assume a scaling behavior of the 
form Lc ^ {F — Fc)~^ where 77 is the associated expo- 
nent. Starting from a flat interface at t = with a fixed 
value of the pushing force F, the horizontal correlation 
length, lc, grows with time whereas the cluster length Lc, 
remains at a fixed value. At some time, that we define 
as the crossover time tc, we have lc{tc) = Lc, and for 
times t > tc the evolution is described by the EW model. 
From its definition, tc has the following scaling relation, 
tc-^ {F - Fc)-", where v = 2ry. 

Regarding the short time behavior for which lc < Lc, 
there is a regime with scaling properties in which the 
exponents and the scaling relations could be defined and 
they would be different from the ones associated with the 
EW model. To determine these quantities, we only need 
to realize that W{t — t') depends on the relevant scale in 
this regime, Lc. However, due to the scaling properties, 
this is equivalent to assuming a dependence on L. In this 
way, we make a very general assumption about the form 
that W{t - t') scales with L and t, W{t - t') - L-^a"* 
where Q, and 7 are two independent exponents. Now, 
substituting this expression of W{t — t') in the effective 
model, we obtain a scaling dependence of cr in L even 
before the saturation time is approached. Hence in this 
case the interface width scales with time and size before 
it saturates (i < tg) as cr ~ t^ instead of the usual 
scaling given bycr - tl^. The expressions for the exponent 
(3, the new exponent a' and the exponent z associated 
with lc as lc ''^ t^ are now given in terms of the parameters 
il and 7 by: 

«' = -77^ P^J^T^ ^ = 2 (7) 
2 + 7 2(2 + 7) 

From Eq.(^, we can discuss the expression of the a 
exponents that one gets for small system sizes for which 
the saturation happens before the effects of the pinning 
clusters have disappeared. The roughening exponent a 
is calculated by taking t — ts L^ in the expression for 
a{t). The result is: 

■i-d n 

a = 7: 7: (8) 

2 + 7 2-f7 ^ ' 
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A model for temporal correlation W(t — t').- In the 
previous paragraph, we characterized the two different 
scahng regimes. Now, we would Hke to propose a simple 
model for W{u, t — t') which match these two regimes. 
Based on physical arguments, one could assume that 
W{u, t — t') is a regular function peaked on the mean 
height u — h{t — t') with a width a. As a simple assump- 
tion, we consider a lorentzian function. In the short time 
behavior it has to have a scaling dependence in L and 
a. This could simply be done by means of the following 
model 

W{u,t)^ " U^^-^ + l) ' (9) 

where 7 and 17 are the parameters that characterized the 
model. The two different scaling regimes are defined by 
L^^ > l?a-^-^ and <C h^cj-^-f. From Eq. (§) we 
obtain that for long time, t ^ tc the scaling behavior of 
W{0, t) is W[0, \t) ~ A-M(t), as A ^ 00, and W{t - t') 
reduces to a delta function, as expected. By substituting 
W{t — t') in the effective model one could get the expo- 
nents in the different scaling regimes. The advantage in 
assuming an expression for W{t — t') over the complete 
temporal regime is that one can determine the crossover 
time, tc, that can be defined hy = h {tc)(T-^-^[t^). 
By substituting h ^ t and the scaling behavior of a one 
obtains 

Numerical results and comparison.- At this point, all 
our exponents and scaling relations associated with the 
first scaling regime depend on the parameters Q and 7. 
To corroborate the presence of both scaling regimes and 
the values of these parameters, we have performed nu- 
merical simulations of Eq.(^ for d = 1 with the same 
parameters as in Ref. In our numerical results, the 
pinning transition occurs at Fc ~ 0.13. We observe a first 
scaling regime with /3 ^ 3/4 in accordance with Parisi re- 
sults. This result implies that our parameter 7 is zero in 
this regime. After a crossover t > tc, (3 ^ 1/4: like the 
EW model. In Fig. 1 we present our numerical results for 
Ic versus time in logaritmic scales. We observe the same 
diffusive behavior for the complete temporal regime in 
accordance with our model. 

Regarding the roughening exponent a, our results are 
in agreement with those obtained by different methods, 
a = 0.7 - 0.9 for ts < t^ ^^j^ and a = 0.49 for t^ > 

Furthermore, in the first regime, the results /3 = 3/4 
and a = 0.8 [Q are compatible in our analytical model 
by taking for our free parameters the values 7 = and 
r2 = 1.4. Finally, in Fig. 2 we present our numerical re- 
sults for the crossover time tc versus F — Fc- We obtain 
an exponent v ^ 1.43, which in terms of L is equivalent 
to tc ~ L^~'^ = L°-^^. This numerical result is in good 



agreement with our analytical results obtained by substi- 
tuting in Eq.(|l^ the values of 17, /3 and 7, i.e. tc ^ L^-^^. 

In conclusion, our analytical results predict a growth 
for the width in agreement with different simulations and 
experiments with the appropriate choice of two free pa- 
rameters . A scaling relation determines the crossover 
time exponent. In our model, the correlation length 
shows a diffusive behavior and we believe it has to be ana- 
lyzed in other models, like directed percolation (DP) and 
the random field Ising model, in order to clarify the diver- 
sity of growth behaviors and hence universality classes. 
For instance, in DP |^ anisotropy generates KPZ expo- 
nents after the crossover instead of those of EW. 
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overhangs because the partial differential equation has a 
single-valued solution. 

FIG. 1. Horizontal correlation length versus time for two 
values of F, close F = 0.15 and far F = 0.4 from the pinning 
transition Fc = 0.13. The slope of the adjusted curve is 1/2. 

FIG. 2. Crossover time versus F — Fc in log-log scale. The 
adjusted curve gives an exponent u = 1.43. 
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